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SUMMARY 


A theory for very high supersonic flow is applied to determine a 
first approximation to the pressure distribution, lift, and drag due 
to lift for a slender cylindrical body of revolution and a wing-body 
combination. The method is applicable only if the altit ude of flight 
is low enough for the gas dynamics type of flow to exist, if the angle 
of attack is large, and if the Mach number is extremely large. The 
•theory used was derived on the basis of a ratio of specific heats equal 
to unity. 

The method indicates that the flow effectively separates from the 
body before reaching the widest part of the cross section. The boundary 
of this region is determined. Certain load concentrations are found 
near the wing-body Junction of the configuration Investigated. 


INTRODUCTION 


Long-range rocket-propelled missiles sometimes reach very high 
angles of attack at high supersonic speeds as they re-enter the dense 
lower atmosphere. These high angles of attack may be reached accidentally 
because of the lack of control at higher altitudes or because of the 
lift required to level out the flight path and to zoom the rocket upward 
again. 

The most efficient shape for obtaining lift at very high supersonic 
speeds is a body with a flat lower surface as pointed out in reference 1. 
The pressure distribution for such a body is the same as that of a flat 
plate and can be calculated by the method of reference 2. A slight 
extension of this method is necessary, however, before the lift of the 
more common type of missile (one with circ ular cross section) c an be 
determined . 
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A method is presented in reference 2 for approxl.mfl.tely predicting 
the drag of three-dimensional shapes and the lift and drag of two- 
dimensional shapes at very hi$i supersonic speeds. An extension of the 
previous work to include the lift and drag of slender cy lin drical 
bodies at high~ angles of attack is given herein. More specifically, 
the present paper is confined to a consideration of the forces on that 
part of a circular cylinder which is not influenced by the nose of 
the body and to a discussion of the spanwise load distribution of a wing- 
body combination at high angles of attack. Determination of the 
pressures on the nose and the effect of nose pressures on body pressures 
are left for future research. 


S2MB0IS 


Cjj drag coefficient due to lift based on plan area of cylinder 

Cj^ lift coefficient based on plan area of cylinder 
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normal-force coefficient based on plan area of cylinder 
constants in equation for the shape of the shock wave 

drag 

diameter 

lift 

length of section of body under consideration 
momentum 

pressure coefficient; any point 
radius 

free -stream velocity 
coordinate axes 
angle of attack 

angle between surface, or shock wave, and free-stream 
direction 
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ratio of specific heats 
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P 

Subscripts : 
1 
2 
h 
i 
N 
n 

V 


free -stream density 

point where high-density region separates from body 
point -where high-density region hits wing 
horizontal 

end point of integration 
normal 

based on normal velocity component 
vertical 


METHOD OF ANALYSIS 


Figure 1 is a schematic representation of a missile, with the 
hatched area being that suitable for investigation by the present method. 


Mach lines 



Figure 1.- Missile with region of application of 
present method (hatched area) • 
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A method for determining the pressures over several aerodynamic shapes 
at very high supersonic speeds in the gas -dynamics flow regime is given 
in reference 2. The method is based on the fact that the pressure 
coefficient at any point on a body approaches a limiting value as the 
Mach number is increased. As the speed increases, the shock wave wraps 
more tightly about the body, and the region of disturbed flow around the 
body becomes relatively small. Keference 2 showed that the flew pattern 
and, hence, velocities, centrifugal forces, and surface pressures are 
not affected greatly by the value chosen for the ratio of specific 
heats 7 . Since the selection of 7 = 1.0 simplifies the problem 
considerably without modifying the results very much, this value of 
7 is used in the present analysis. 

Figure 1 clearly shows the similarity between a cylinder at high 
angle of attack and a wing swept back outside the Mach cone. As on a 
sweptback wing, the pressures on the cylinder are determined only by 
the component of the free -Btr earn velocity normal to the axis of the 
cylinder. The combination of free -stream Mach number and angle of attack 
must be such that the component % is much greater than the local 

speed of sound in order for the calculations to be valid. For a given 
forward Mach number, therefore, the equations are more accurate at high, 
angles of attack- The problem has thus been reduced to the investigation 
of the pressures over a circular cylinder normal to the stream of air 
moving at Tjf. 

A schematic representation of the general flow pattern based on 
the results of reference 2 is given in figure 2. The thickness of the 
high-pressure region has been exaggerated for clarity. Also, the flow 
is shown to separate from the surface at some point which is to be 
calculated subsequently in this section- 
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Shock vave 



la figure 3 . ai’e shown the angle and points that are used in the 
calculations . 



Figure 3 .- Geometry of the analysis 
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The pressure coefficient P at any point (xi,yj_) is given in reference 2 
as 


ap ± py i 

P_. = 2 sin%4 + 2 sin P^ - — cos P dy (l) 

1 ^i Jo 


For the circular cylinder, ^ = cos p and = - — ' H .j ~p • Then, the 
pressure coefficient for the circular cy lind er at any point is 


p± . 2(1 


cos 2 P^) 


- 2-siUi F 7i cos f 

r sin P ± J 0 


< 3 y 




= 2 - 



( 2 ) 


A sketch of the pressure distribution over the cylinder is given 
as figure 4; a more detailed plot of the pressure distribution is given 
in figure 5- 



Figure 


k.- Pressure distribution over the cylinder- 
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The limi ting negative pressure coefficient at very high supersonic speeds 
is quite close to zero, and. for practical purposes negative pressure 
coefficients can he assumed, to he unattainable at these speeds. The 
centrifugal force due to the curvature of the flow is equal to the 
difference in pressure between the shock wave and the surface of the 
cylinder. When the calculated surface pressure coefficient reaches zero, 
the shock wave mas t diverge from the body inasmuch as the pressure drop 
across the disturbed air layer is no longer sufficient to produce the 
flow curvature req uir ed for following the surface. The separation point 
(xi,yi) where the shock wave can no longer follow the surface curvature 
is found by setting the pressure coefficient equal to zero; that is. 



therefore. 


yi 

r 


~ = o .8165 


(3) 


Since for a body with a given diameter the lift and drag increase 
directly with the length of the body, the force coefficients are based 
on the plan area of the cylinder - that is, the product of length and 
diameter. The normal-force coefficient based on the dynamic pressure 
corresponding to the normal velocity component Vpj is 




'll 


21 


p i fly 


Zd 




= 1.089 


(4) 


Based on free -stream, dynamic pressure, the normal-force coefficient is 


°N = % I f j = i* 08 ? 


(5) 
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Then, the lift coefficient is 


= I.089 sin^a, cos a 


( 6 ) 


and the drag coefficient due to lift is 


Cp = I.O89 sin^a 

The lift-drag ratio, therefore, is 


( 7 ) 


L 

D 


cot a 


( 8 ) 


Equations (6), (7), and (8) are plotted in figures 6, 7 , and 8, respectively . 
Maximum, lift is obtained at an angle of attack of about 55 ° (fig* 6); 
the corresponding drag is very high (fig. 7 ) • Because of the neglect 
of viscous drag, very high lift-drag ratios are indicated for low angles 
of attack (fig. 8) . The method should not, however, be applied at low 
angles. 

The shape of the shock wave that wraps about the cylindrical body 
can be estimated from momentum considerations. The "zero -pressure stream- 
line method" of reference 2 assumes that the very -high-density region 
behind the shock wave is very narrow and is bounded on one side by 
shock pressure and on the other side by negligible pressure. Since the 
shock wave is one of the boundaries of the, very -narrow, high-density region, 
the direction of the shock wave at any point is the same as the direction 
of the total momehtum in the high-density region at that same pointy 
therefore, the slope of the shock wave at any point (xj; ,y^) is 



The vertical momentum at the point where the flow effectively 

separates from the body is given by the following equation: 


= P"Vjj sin Pi 


^1 

Jo 


cos P dy 
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The horizontal momentum at this same point is given as 




Mhp = pVjj cos Pp 


cos P dy 


UO 


If no further loss in momentum is assumed to take place after the 
high-density flov separates from the "body, the vertical momentum remains 
constant and equal to the value at the point (xp,yp) . The horizontal 

momentum, however, at some point beyond the point of separation 

is the sum of the horizontal momentum at the point (xp,yp) and the 

horizontal momentum added from the free stream between the points 
(xp,yp) and (x^yp) • Thus, 


*i 


Mhp = %j_ + j dy 


2 / P 7 1 

= PVjf cos Pp cos P dy - yp + yp 

do 


Substituting the foregoing relations for and in 

equation (9) and letting 

' P7l 

Cp = 2 ( cos Pp cos P dy - yp 

Jo 


and 


17j_ 

Cg = -2 sin Pp cos P dy 

do 
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give 


= ~ C 2 

c i + 2y ± 


Integrating yields the following relation: 


7i 2 + C iyi + = C 3 


( 10 ) 


The constant may he obtained hy evaluating the left side of the 

equation at the point (xq, yq) * The equation of the shock wave- is, thus, 
parabolic, and the constants Involved vary with the shape of the body. 

For the case considered herein (cylinder with r = 1.0 ) the 
constants xn&y be evaluated from the following relations (see fig. 3) ; 


sin P = 1 - x 


cos P = y 



7q 

cos P dy = 


Slnce yq was found to.be O .8165 (equation (3)), then 


tq = r - ^r 2 - yq 2 = OA 23 I 


Evaluating the constants and substituting them in equation (10) 
yield the following equation for the shock wave after it leaves the 
cylinder: 


y ± 2 - 1.0897! " 0 . 3852 ^ + O .385 = 0 


( 11 ) 


Equation (ll) is plotted in figure 9* 



MCA 03? No . 17^0 


11 


SPANWISE LOAD DISTRIBUTION 


The method used herein for obtaining the spanwi.se load distribution 
is limited to wing sections with flat lower surfaces. The span-wise load 
distribution can then be calculated in three parts. 

Load distribution over cylindrical body .- The load distribution over 
the cylindrical body can be obtained from equation (2) or from figure 5* 

Load distribution along wing outboard from point where high-density 
air hits wing .- Since a complete loss in momentum fin the direction normal 
to the body axis) occurs on the lower surfaces of the wings just as on 
the leading edge of the cylindrical body, the pressure coefficient is, 
therefore, the same at both places. This pressure coefficient is the 
peak pressure coefficient of the cylindrical body, and the value is 2 .0 . 

Concentrated load in high-density flow that separates from cylinder . - 
The pressure near the wing-body intersection follows from elementary 
considerations. Since the normal-force coefficient (l.089)-of the body 
based on the dynamic pressure for the velocity normal to the axis of the 
body is less than that due to a complete loss in momentum. (2.0), 
appreciable momentum remains in the air stream. This momentum is 
concentrated into a very narrow span-wise distance (zero for the case of . 
7 - 1.0 because the density goes to infinity) and, finally, is completely 
lost at the wing surface (near the wing root) . The momentum of the flow 
normal to the axis of the body is a l 1 lost (part at the body and part 
near the wing-body Juncture); therefore, the total lift for a unit 
length of wing and body is the same as that for a continuous wing. 

Inasmuch as the lift coefficient over the body is about one -half that 
over a wing, the lift coefficient at the wing root must then be higher 
than that of a wing alone; thus, very high peak pressures exist near 
the wing root. The theory indicates that for the limiting case of 
extremely high Mach numbers the pressure coefficient is infinity for 
an infinitesimal distance along the span (point Po in fig. 10). Point Pg 
can be found from the equation of the shock wave (equation (ll)). for 
finite Mach numbers the peak pressures are expected to be finite and 
the pressure distribution is expected to have less abrupt changes . 

The part of the body ahead of or behind the section influencing the 
wing develops lift in the manner shown in figures 5 and 6. 


CONCLUDING- REMARKS 


An approximate theory for very high supersonic speeds is applied 
to determine the lift and drag due to lift for a cylindrical body and a 
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■wing-body combination at high angles of attack in a gas dynamics type 
of flov- For a body with a given diameter the lift is directly proportional 
to the body length. The lift of the circular body is approximately one- 
half that of a wing of the same plan area. The wing roots have extremely 
high pressures because they convert the flow around the local body 
sections into additional lift. The wing and the corresponding part of 
the body develop a combined, lift equal to that of a continuous wing. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., August 3 , 1948 
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Pressure dreg coefficient due to 
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Figure 7*“ Pressure drag ooeffioient for seotion of oylinder investigated. 
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Ratio of lift to pressure drag due to 



figure <5.“ Ratio of lift to pressure drag due to lift for seotlon of oyllnder Investigated 
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figure 9»” Oaloulated shocfc-wave shape. 
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Figure 10.— Bpanwlae load distribution for wing-body oombinatlon 



